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Three-body recombination is a prime example of the fundamental interaction between three particles and it is of importance to the 
physics of ultracold gases. Due to the complexity of this process it has resisted a comprehensive theoretical description. Experimental 
investigations have mainly focussed on the observation of corresponding loss rates without revealing information on the reaction 
products. Here, we provide the first general experimental study on the population distribution of molecular quantum states after 
three-body recombination. We utilize a novel detection scheme which combines photoionization of the molecules with subsequent ion 
trapping. By analyzing the ionization spectrum, we identify the population of energy levels with binding energies up to h x 750 GHz. 
We find a broad population of electronic and nuclear spin states and determine a range of populated vibrational and rotational states. 
The method presented here can be expanded to provide a full survey of the products of the recombination process. This may be pivotal 
in developing an in-depth model of the underlying three-body dynamics. 



While cold collisions of two atoms are understood to an excellent 
degree, the addition of a third collision partner drastically complicates 
the interaction dynamics. In the context of Bose-Einstein condensation 
in atomic gases, three-body recombination plays a crucial role^' and 
it constitutes a current frontier of few-body physics^'. However, the 
investigations focussed mainly on the atom loss rates established by 
the recombination events. Discussions of the final states populated in 
the recombination process were restricted to the special case of large 
two-body scattering lengths^ a nd cu lminated in the prediction and 
observation of Efimov resonance ;! 10 * 12 ! In general, recombinatio n pro- 
cesses are of fundamental interest in various physical systems^ 13 * 14 ! 
The control and tunability of ultracold atomic systems provides an ex- 
perimental testbed for a detailed understanding of the nature of these 
processes. 

Here, we demonstrate the probing of molecules with binding ener- 
gies up to h x 750 GHz (where h is Planck's constant) generated via 
three-body recombination of ultracold thermal 87 Rb atoms. We pro- 
duce the atomic sample in an optical dipole trap located within a linear 
Paul trap. The recombination and detection process is illustrated in 
Fig. la-d. Following a recombination event, the created Rb2 molecule 
can undergo resonance-enhanced multi-photon ionization (REMPI) by 
absorbing photons from the dipole trap laser at a wavelength around 
1064.5nm. The ion is then captured in the Paul trap and detected essen- 
tially background-free with very high sensitivity on the single particle 
level. Fig. le shows a simplified scheme of the Rb2 and Rb^ poten- 
tial energy curves. From weakly-bound molecular states three photons 
suffice to reach the molecular ionization threshold. An additional pho- 
ton may dissociate the molecular ion. By scanning the frequency of the 
dipole trap laser by more than 60 GHz we obtained a high resolution 
spectrum featuring more than 100 resonance peaks. This dense and 
complex spectrum contains information which vibrational, rotational 
and hyperfine levels of the Rb2 molecule are populated. We present an 
analysis of these data and make a first assignment of the most promi- 
nent resonances. This assignment indicates that in the recombination 
events a broad range of levels is populated in terms of vibrational, ro- 
tational, electronic and nuclear spin quantum numbers. 

Our experimental scheme to detect cold molecules makes use of 
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Figure 1 | Illustration of recombination and ionization in the 
atom-ion trap, a, A three-body collision in the ultracold gas of 87 Rb 
atoms leads to a recombination event in which, b, a Rb2 molecule is formed 
with high kinetic energy. C, While the atom is lost from the trap, the 
molecule can be photoionized in a REMPI process and trapped in the Paul 
trap, d, The relative positions of the atom and ion trap centers are shifted 
by about 300 fim to avoid atom-ion collisions, e, Potential energy curves of 
the Rb2 and RbJ molecule adapted from refs l 15 i 16 l The curves A, b, c, are 
A Ej", b 3 n u , c 3 E+. The internuclear distance is given in units of Bohr 
radii ao. A REMPI path with three photons is shown. It can create RbJ 
ions in vibrational states up to v ss 17. 
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Figure 2 | Dependence of the ion production rate r ioll on atomic 
density. Fi on is normalized to a cloud atom number of 10 6 atoms. The data 
are well described by a quadratic fit (solid green line). They were taken at a 
constant dipole trap laser intensity II = 36 kW/cm 2 . 

the generally excellent detection efficiencies attainable for trapped 
ions. It is related to proven techniques where photoassociated cold 
molecules in magneto-optic al traps were photoionized from the sin- 
glet and triplet ground state d 1 ^ 17 -^ (see also ref^). Our method is 
novel as it introduces the usage of a hybrid atom-ion trap which signif- 
icantly improves the detection sensitivity. We perform the following 
experimental sequence. A thermal atomic sample typically containing 
iVat ~ 5 x 10 5 spin-polarized 87 Rb atoms in the \F = l,mp = — 1) 
hyperfine state is prepared in a crossed optical dipole trap at a magnetic 
field of about 5 G. The trap is positioned onto the nodal line of the ra- 
diofrequency field of a linear Paul trap. Along the axis of the Paul trap 
the centers of the atom and ion trap are separated by about 300 um to 
avoid unwanted atom-ion collisions (Fig. Id). At atomic temperatures 
of about 700nK and peak densities no ~ 5 x 10 13 cm -3 the total three- 
body recombination rate in the gas is r rec = I/371 2 A r al /3 5//2 ~ 10kHz. 
Here, the three-body loss rate coefficient L3 was taken from ref^. At 
the rate r rec , pairs of Rb2 molecules and Rb atoms are formed as final 
products of the reactions. Both atom and molecule would generally be 
lost from the trap due to the comparatively large kinetic energy they 
gain in the recombination event (in our case typically on the order of a 
few K xfce where fee is the Boltzmann constant). The molecule, how- 
ever, can be state- selectively ionized in a REMPI process driven by 
the dipole trap laser. These molecular ions remain trapped in the deep 
Paul trap and are detected with single particle sensitivity (see Meth- 
ods). In each experimental run, we hold the atomic sample for a time 
r ss 10 s. After this time we measure the number of produced ions 
in the trap from which we derive (after averaging over tens of runs) 
the ion production rate Tion normalized to a cloud atom number of 10 6 
atoms. 

As a consistency check of our assumption that Rb2 molecules 
are ionized in the REMPI process, we verify the production of RbJ 
molecules. For this, we perform ion mass spectrometry in the Paul 
trap (see Methods). We detect primarily molecular RbJ ions, a good 
fraction of atomic Rb + ions but no RbJ ions. Our experiments show 
that Rb + ions are produced in light-assisted collisions of Rb^ ions 
with Rb atoms on timescales below a few ms. Details of this dissocia- 
tion mechanism are currently under investigation and will be discussed 
elsewhere. 

In the following, we analyze the dependence of the ion production 
rate on atomic density as well as on laser frequency and intensity. This 
will ultimately help understand the three-body recombination process. 



By investigating the dependence of the ion production rate on 
atomic density, we verify that three atoms are originally involved in 
the reaction. For this measurement, the density is mainly adjusted by 
varying the cloud atom number while keeping the light intensity of the 
dipole trap constant. We find the obtained data (Fig. 2) to be in good 
agreement with a fit curve assuming a simple quadratic dependence 
Tio,, oc tiq. Thus, a three-body process marks the starting point of the 
ionization process. 

Next, we investigate the dependence of the ion production rate on 
the wavelength of the narrow-linewidth dipole trap laser (see Meth- 
ods). We scan the wavelength over a range of about 0.3 nm around 
1064.5 nm, corresponding to a frequency range of about 60 GHz. Typ- 
ical frequency step sizes are 50 MHz or 100 MHz. We obtain a rich 
spectrum of resonance lines which is shown in Fig. 3a. The quantity 
Fion denotes the ion production rate normalized to the atom number of 
the cloud and to the square of the atomic peak density. 

We find strongly varying resonance strengths and at first sight 
fairly irregular frequency spacings. In the following we will argue 
that most resonance lines can be attributed to respective well-defined 
molecular levels (resolving vibrational, rotational and often even hy- 
perfine structure) that have been populated in the recombination pro- 
cess. These levels are located in the triplet or singlet ground state, 
a 3 Sj and X 1 ^ , respectively. The relatively dense distribution of 
these lines reflects that a fairly broad range of states is populated. 
A direct assignment of the observed resonances is challenging, as it 
hinges on the precise knowledge of the level structure of all the rel- 
evant ground and excited states. In the following we will access and 
understand the data step by step. 

One feature of the spectrum that catches the eye is the narrow 
linewidth of many lines. For example, Fig. 3b shows a resonance of 
which the substructures have typical half-widths Az/ r w 50MHz. This 
allows us to estimate the maximal binding energy of the molecules 
involved. Since the velocity of the colliding ultracold atoms is ex- 
tremely low, the kinetics of the recombination products is dominated 
by the released molecular binding energy E\>. Due to energy and mo- 
mentum conservation the molecules will be expelled from the reac- 
tion with a molecular velocity «Rb2 = \/2Eb/ ('immi) where rriKbi 
is the molecular mass. The molecular resonance frequency Vo will 
then be Doppler-broadened following a distribution with half-widths 
Ai/D = y/^vottssa/^c Here, c is the speed of light. By comparing 
Az^d to the observed values of Av r we estimate a maximal binding 
energy of i?b,max 2.5 THz. This simple analysis overestimates 

the value £b,max since it neglects the natural linewidth of the transition 
and possible saturation broadening. Still, it already strongly constrains 
the possible populated molecular levels that are observed in our exper- 
iment. 

Next, we investigate the dependence of the ion production rate 
on laser intensity II- In our experimental setup, this measurement is 
rather involved because the laser driving the REMPI process also con- 
fines the atomic cloud. Thus, simply changing only the laser intensity 
would undesirably also change the density no of the atoms. To pre- 
vent this from happening we keep no constant (no ~ 5 x 10 cm ) 
by adjusting the atom number and temperature appropriately. Due 
to these experimental complications we can only vary II roughly by 
a factor of 2 (Fig. 4a). We set the laser frequency to the value of 
v l = = 281610 GHz, on the tail of a large resonance (see Fig. 3). 
The atomic temperatures in this measurement range between 500 nK 
and 1.1 uK, well above the critical temperatures for Bose-Einstein con- 
densation. The atomic densities can therefore be described using a 
Maxwell-Boltzmann distribution. Assuming a simple power-law de- 
pendence of the form Ti on oc I" we obtain the best fit using an ex- 




Figure 3 | REMPI spectrum, a, A scan of the dipole trap laser frequency over more than 60 GHz around an offset frequency = 281.610 THz 
shows a multitude of resonance lines. Each data point is the result of 30 to 60 repetitions of the experiment with ion detection on the single particle level. 
The total spectrum was obtained over a time span of 2 months. Checks of the long-term consistency of resonance positions and strengths were performed. 
Spectral regions dominated by transitions to c 3 E+ are indicated by the shaded areas in dark and light blue (0~ and l g component, respectively), b, 
High-resolution scan of the strong resonance at — i/^ « 0.5 GHz. c, Central spectral region with assigned P/R branches of the transition X 1 E+ (v = 
115) — >• j4 1 E+(i/ = 68). The corresponding quantum numbers I and J' are given). P(J) marks the transition J — ► J+l, R(J) the transition J — > J — l. 
The region where also transitions to 6 3 n u appear is marked by a green horizontal bar. 



ponent a = 1.5(1) (solid green line in Fig. 4a). This fit is between a 
linear and a quadratic intensity dependence (dashed red and blue lines, 
respectively). Thus, at least two of the three transitions composing the 
ionization process are partially saturated at the typical intensities used. 

In order to better circumvent possible density variations of the atomic 
cloud induced by changes in laser intensity, we employ a further 
method which enables us to vary the intensity with negligible effects 
on the atomic sample. We achieve this by keeping the time-averaged 
intensity (7l) constant and comparing the ion production rates in a 
continuous dipole trap with a "chopped" dipole trap in which the in- 
tensity is rapidly switched between and 2/l. In both cases the trap 
is operated at an intensity (Jl) ~ 15 kW/cm 2 . In the "chopped" con- 
figuration the intensity is switched at a frequency of 100 kHz so that 
the atoms are exposed to the light for 5 /is followed by 5 /is without 
light. It should be noted that molecules formed in the "dark" period 
with sufficiently high kinetic energies may leave the central trapping 
region before the laser light is switched back on. They are then lost for 
our REMPI detection. Taking into account the molecular velocity and 
the transverse extensions of the laser beams we can estimate that this 
potential loss mechanism leads to errors of less than 30%, even at the 
highest binding energies relevant to this work (Eb ~ h x 750GHz, see 
below). We did not observe evidence of such losses experimentally. 
Investigations were made by changing the chopping frequency. We 
define R as ratio of the ion production rates in the "chopped" and the 
continuous trap configuration. Fig. 4b shows the results of these mea- 
surements for various laser frequencies ul- We find a value R ~ 1.5 
for off-resonant frequency settings — < 0.4 GHz, in good agree- 
ment with the result presented in Fig. 4a. When scanning the laser onto 
resonance at ul — ~ 0.45 GHz (see also Fig. 3b) we obtain R w 1, 
This result indicates a linear intensity dependence of the REMPI pro- 



cess in the resonant case, which is explained by the saturation of two 
of the three molecular transitions involved. It is known that transitions 
into the ionization continuum (photon III, see Fig. le) will not satu- 
rate under the present experimental conditions. This means that the 
excitation pathway via photon I and II must be saturated and there- 
fore on resonance. However, an inspection of the level structure shows 
that photon I can only resonantly reach three vibrational levels which 
correlate to the 5s + 5p asymptote. This is largely due to its long 
wavelength of 1064.5 nm. For example, within the c E+ state only 
the lowest («' = 0) vibrational level can be excited, namely from level 
v = 26 of a 3 E+ which has a binding energy of Eb ~ h x 516 GHz. 
The other two excited vibrational levels are the v' = 68 of A ll S u ,0i 
and v = 72 of b 3 H u , 0^. The dipole trap laser couples these levels 
resonantly to v = 115 of X 1 ^ (binding energy ~ 161 GHz x h) and 
v = 109 of A^Ejj" (binding energy « 754 GHz x h), respectively (see 
Fig. 5 for details). Coupling to other nearby a or X ground state lev- 
els can be neglected due to weak transition strengths, resulting from 
small Franck-Condon factors or dipole selection rules. We note that 
there is a strong mixing of A 1 ^ and b s Tl u due to spin-orbit coupling. 
The transitions that we consider exhibit large enough Franck-Condon 
factors (typ. 10~ 2 ...10 -3 ) so that at laser powers of ~ 10 4 W/cm 2 
resonant transitions can be well saturated. 

From recent spectroscopic studies^^ and additional measure- 
ments in our lab 25 the level structure of all relevant levels of the 
X, a, A, b, and c states is well known. The absolute precision of most 
of the level energies is far better than 1 GHz for low rotational quan- 
tum numbers J. Spin-orbit and effective spin-spin coupling in the A, 
b, and c states lead to Hund's case c coupling where the relevant levels 
of states A and 6 have U symmetry while the levels of state c 3 Ej" 
are grouped into 0,7 and 1 9 components. The level structure of the 0^ 
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Figure 4 | Dependence of the ion production rate on the intensity of the dipole trap laser, a, Assuming a power-law dependence f ion oc I", 
the best fit to the data is achieved for a « 1.5 (solid green line). Linear and quadratic fits are also given (blue and red dashed lines, respectively), b, 
Measurement of the intensity dependence using a "chopped" dipole trap. The ratio R w 1 on resonance indicates saturation of both transitions I and II. 



states is quite simple as it is dominated by rotational splittings. Typical 
rotational constants for the electronically excited states are on the or- 
der of 400 MHz, for the weakly bound X and a states they are around 
100-150 MHz. 

In the experimental data (Fig. 3a) the central region from vl — 
l/£ = —6 to 7 GHz is marked by several prominent resonances that 
are significantly stronger than those observed throughout the rest of 
the spectrum. These resonance peaks can be explained by transitions 
from the X ground state to A and b states. The prominence of these 
singlet transitions is explained by the degeneracy of levels due to small 
hyperfine splittings. Indeed, by analyzing these strong resonances 
with regard to line splittings and intensities it was possible to consis- 
tently assign rotational ladders for total nuclear spin quantum numbers 
I = 1, 2, 3 for the transition A^E+O = 115) ^E+O' = 68). 
For the weakly bound X states, the / quantum numbers are approxi- 
mately good. For 7 = 1,3 the total parity of the molecule is negative, 
while for I = 0, 2 the parity is positive. The starting point of the 
rotational ladder is fixed by the already mentioned additional spectro- 
scopic measurements 2 ^! In Fig. 3c the P branch and the R branch of 
the rotational ladders are indicated. These lines can be grouped into 
pairs sharing the same J' of the excited state and I quantum number. 
At frequencies — z^£ > 2 GHz additional strong lines appear that 
we attribute to the X l H+(v = 109) -> b 3 U u (v' = 72) transition. 
We can now roughly estimate the range of molecular rotation J of 
the populated levels in the ground state. The strong isolated lines that 
we have assigned to the A X E^ — > transition are all contained 

within a relatively small spectral region 

(K - v °\ < 6 GHz) and 
are explained by rotational quantum numbers J < 7. Population of 
higher rotational quantum numbers would result in a continuation of 
the strong resonance lines stretching to transition frequencies beyond 
vl — z^£ = 10 GHz, which we do not observe. Similarly, if only ro- 
tational quantum numbers J < 5 were populated, a spectrum would 
result which does not have enough lines to explain the data. Thus, we 
can roughly set the limits on the molecular rotation to J < 7. These 
results are consistent with our observations of the spread of the transi- 
tions X — > b and a — > c (see Fig. 5). 

From our measurements we can also make an estimate about the 
number of molecular vibrational levels populated in the recombination 
events. For this we integrate over the spectral region dominated by the 
transition X x E+(?; = 115) -> A 1 E+(w' = 68) and divide by the 
observed FWHM linewidth (w 100 MHz). The result is then divided 
by a factor of 2 to account for the P and R branch of the transition. For 
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Figure 5 | Overview over relevant molecular transitions. The en- 
ergy E of the energetically lowest levels of each vibrational manifold is 
given with respect to the 5s5s asymptote. Colored thick arrows represent 
molecular transitions relevant to the spectrum of Fig. 3. The expected rel- 
ative strengths of these transitions are also given. Here we only consider 
Franck-Condon factors and mixing of singlet and triplet, while electronic 
transition moments are ignored. Grey arrows mark transitions that oc- 
cur in the relevant spectral region but are very weak. We identify three 
main molecular transitions for the initial step of the REMPI process. (1) 
Molecules in the v = 26 vibrational level of the a 3 S+ potential are excited 
to the v' = level of the c 3 E+ potential. This level is split into a l g and 
a 0~ component. (2) Excitation from A 1 E+(i> = 115) to J 4 1 E+(^' = 
68). (3) Excitation from A^E+O = 109) to 6 3 n u (V = 72). Transition 
(3) becomes possible through the strong spin-orbit coupling of the A and b 
states. 



the given atom numbers and atomic densities in our cloud, this yields 
an estimate for the global ion production rate from this vibrational 
level. Assuming that the two first transitions of the REMPI process 
are saturated and using a typical ionization cross section for the third 
transition to be ~10 MBanP^I we obtain an ionization probability 
of around 10 % for molecules formed in one of the levels susceptible 
to the REMPI process. From this we can estimate the population rate 
for the X 1 'E£(v = 115) level due to three-body recombination. It is 
roughly on the order of 0.2 kHz. We can compare this to the much 
larger total recombination rate of r rec = 10 kHz in the atomic gas. 
Consequently, our estimate suggests that on the order of 50 vibrational 
levels, a significant fraction of all the 169 existing levels of the X and 
a states, may be populated in the three-body recombination. This esti- 
mate seems reasonable considering that we have observed molecule 
production down to X 1 E^(u = 109), which is located 38 vibra- 
tional levels (counting both singlet and triplet states) below the atomic 
asymptote. 

In conclusion, we have carried out a first, detailed experimental study 
of the molecular reaction products after three-body recombination of 
ultracold Rb atoms. We use a high-power, narrow-linewidth laser to 
state-selectively ionize the produced molecules in a REMPI process. 
Subsequently, these ions are trapped in an ion trap and detected with 
very high sensitivity and negligible background. An analysis of the 
ionization spectrum allows us to identify population of several vibra- 
tional quantum levels indicating that the recombination events result 
in a fairly broad and uniform population distribution. We conjecture 
that dozens of vibrational levels are populated in total. Molecules are 
produced both in X 1 ^ as well as a 3 £+, with negative and positive 
total parity, various total nuclear spins and rotational quantum numbers 
reaching J < 7. We have made a first assignment of several strong in- 
dividual spectral lines. Confirming and expanding this assignment to 
other transitions will require further experimental and theoretical ef- 
forts to acquire a more detailed knowledge of the involved molecular 
levels. Such studies may finally pave the way to a detailed understand- 
ing of three-body recombination in the future. 

Reaching beyond the scope of three-body recombination, the great 
sensitivity of our detection scheme has enabled us to state-selectively 
probe single molecules that are produced at rates of only a few Hz. 
We thereby demonstrate a novel scheme for precision molecular spec- 
troscopy in extremely dilute ensembles. 

METHODS 

Dipole trap and REMPI configuration. The crossed dipole trap is com- 
posed of a horizontal and a vertical beam focussed to beam waists of 
~90 /im and ~150 ^m, respectively. It is positioned onto the nodal line 
of the radiofrequency field of the linear Paul trap with fj,m precision. The 
two trap centers are separated by about 300 /im along the axis of the Paul 
trap (see Fig. Id). In a typical configuration, the trap frequencies of the 
dipole trap are (175, 230, 80) Hz resulting in atom cloud radii of about 
(6, 7, 16) fim. The short-term frequency stability of the dipole trap laser 
source is on the order of 1 kHz and it is stabilized against thermal drifts to 
achieve long-term stability of a few MHz. The two beams of the dipole trap 
are mutually detuned by 160MHz to avoid interference effects in the optical 
trap. Consequently, two frequencies are in principle available to drive the 
REMPI process. However, the intensity of the horizontal beam is 4 times 
larger than the one of the vertical beam and we have not directly observed a 
corresponding doubling of lines. Further details on the atom-ion apparatus 
are given in refl22l 

Paul trap configuration. The linear Paul trap is driven at a radiofrequency 
of 4.17 MHz and an amplitude of about 500 V resulting in radial confine- 
ment with trap frequencies of (a) x ,Bai ^y.Ba) = 2n x (220, 230) kHz for 
a 13S Ba+ ion. Axial confinement is achieved by applying static voltages 



to two endcap electrodes yielding w z-Ba = 2-7T x 40.2 kHz. The trap fre- 
quencies for "dark" RbJ and Rb+ ions produced in the REMPI processes 
are (m Ba /m dark x w^, m^/m^ x w y-Ba , \/m Ba /m dark x w z-Ba ) where 
m Ba and m aar k denote the mass of the Ba+ ion and the dark ion, respec- 
tively. The depth of the Paul trap depends on the ionic mass and exceeds 
2 eV for all ionic species relevant to this work. 

Ion detection methods. We employ two methods to detect Rb J and Rb+ 
ions both of which are not amenable to fluorescence detection. In the first 
of these methods we use a single U'apped and laser-cooled 138 Ba+ ion as 
a probe. By recording its position and trapping frequencies in small ion 
strings with up to 4 ions we detect both the number and the masses of the 
ions following each REMPI process (see als(p2l). The second method is 
based on measuring the number of ions in the Paul trap by immersing them 
into an atom cloud and recording the ion-induced atom loss after a hold 
time of 2 s (see also^). During this detection scheme, we take care to 
suppress further generation of ions by working with small and dilute atomic 
clouds and by detuning the REMPI laser from resonance. Both methods are 
background-free in the sense that no ions are captured on timescales of days 
in the absence of the atom cloud. Further information on both detection 
methods is given in the Supplementary Information. 
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Supplementary Information 

In this Supplementary Information we describe two methods that 
we employ to detect small numbers of Rb^ and Rb + ions in our linear 
Paul trap. 

Method 1 

To implement our first ion detection method allowing mass-sensitive 
detection of "dark" ions we rely on the presence of a single "bright" 
ion in the trap. Information on additional ions can be extracted from 
its fluorescence position. When using this method, our experimental 
procedure begins with the loading of a single 138 Ba + -ion into our 
linear Paul trap. We laser-cool the ion and image its fluorescence 
light onto an electron-multiplying charge-coupled device camera. 
This enables us to determine the position of the trap center to 
better than 100 nm. The ion is confined at radial and axial trapping 
frequencies cj r ,B a ~ 2n x 220 kHz and w ra ,Ba ~ 2-7T x 40.2 kHz and 
typically remains trapped on timescales of days. Next, we prepare an 
ultracold atomic sample in the crossed dipole trap. At typical atomic 
temperatures of about 700 nK the atom cloud has radial and axial 
extensions of about 7 and 15 /im and is thus much smaller than 
the trapping volume of our Paul trap. To avoid atom-ion collisions 
we shift the Ba + -ion by about 300 /im with respect to the atom cloud 
before the atomic sample arrives in the Paul trap. The shifting is 
performed along the axis of the trap by lowering the voltage on one 
of the endcap electrodes. Additionally, we completely extinguish all 
resonant laser light so that the atoms are only subjected to the light 
of the dipole trap. The atomic sample is moved into the center of the 
radial trapping potential of the Paul trap and is typically held at this 
position for a time Thoid ~ 10 s. Despite the axial offset from the 
center of the Paul trap, the atom cloud at this position is fully localized 
within the trapping volume of the Paul trap. After the hold time the 
sample is detected using absorption imaging. Subsequently, the ion 
cooling beams are switched back on for fluorescence detection of the 
Ba + -ion. 

The presence of a second ion in the trap leads to positional shifts of the 
138 Ba + -ion by distances on the order of 10 /im (see Fig.QJ. We make 
use of the mass-dependent trap frequencies of the Paul trap to gain 
information on the ion species trapped. In a two-ion Coulomb crystal 
composed of a Ba + -ion and a dark ion, the axial center-of-mass 
frequency Wax,2ioii shifts with respect to a; ax ,Ba depending on the mass 
of the dark ion mdni 32 . We measure cj ax ,2ion by modulating the 
trap drive at frequencies w mo a and by monitoring the induced axial 
oscillation of the Ba + -ion, visible as a blurring of the fluorescence 
signal. In this way, after each ion trapping event, we identify a 
resonance either at w mo a ~ 27r x 44 kHz or u) m od ~ 2n x 38 kHz corre- 
sponding to md al k = 87u and mdark = 174u, respectively (see tableQJ. 



We have expanded this method for ion strings with up to four ions 
including the Ba + -ion. For this purpose, we perform the following 
step-by-step analysis. 

1. The position x of the Ba + -ion with respect to the trap center is 
detected. If x ^ 0, the value of x allows us to directly determine 
the total number of ions in the string. 

2. If x = we need to distinguish between a single Ba + ion and a 
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Figure 1 | Ion detection using a 138 Ba+-ion. Positional shifts of the 
fluorescence of the Ba+-ion and measurements of the trap oscillation fre- 
quencies allow us to perform mass-sensitive detection of up to three "dark" 
ions in the trap. 

Table 1 | Trap oscillation frequencies of two-ion crystals 



Ion species 


^ax.2k.n/27T [kHz] 


w r /27T [kHz] 


laH B a + and las Ba+ 


40.2 


220.0 


13s Ba+and 87 Rb+ 


44.0 


345.3 


138 Ba+ and 87 Rb^ 


37.7 


170.7 



three-ion string with Ba at its center. This is done by modulating 
the trap drive at Wax3a> thereby only exciting the Ba + ion if no 
further ions are present. 

3. We destructively detect the Rb + ions by modulating the trap drive 
on a 5 kHz wide band around 2 x o; r ,Rb/(27r) = 691 kHz. This 
selectively removes only Rb + ions from the string making use 
of the relatively weak inter-ionic coupling when exciting the ions 
radially. 

4. Steps 1. and 2. are repeated to detect the number of remaining 
ions. 

5. The Rb^ ions are destructively detected via modulation around 

2 x uj r _ Rb2 /(2n) = 341 kHz. 



Method 2 

We have also developed a second ion detection method that does not 
require an ion fluorescence signal. Instead, the trapped ions are de- 
tected via their interaction with an atomic sample. For this purpose, we 
produce a comparatively small atom cloud containing about 1 x 10 5 
atoms at a density of a few 10 12 cm -3 . In addition, we set the fre- 
quency of the dipole trap laser to an off-resonant value so that the pro- 
duction of additional ions during the ion probing procedure becomes 
extremely unlikely. We now fully overlap the ion and atom traps for 
an interaction time r; nt = 2 s. By applying an external electric field 
of several V/m we set the ion excess micromotion energy to values on 
the order of tens of fce xmKl 3 ' l 33 [ Consequently, if ions are present in 
the trap, strong atom losses occur due to elastic atom-ion collisions. 
Fig|2]shows a histogram of the atom numbers of the probe atom sam- 
ples consisting of the outcome of about 1,000 experimental runs. The 
histogram displays several peaks which can be assigned to the discrete 
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Figure 2 | Ion detection method based on ion-induced atom loss. 

We overlap an ultracold atom cloud containing approximately 110,000 
atoms with the center of the Paul trap. After an interaction time t = 2 s 
we detect the ion-induced atom loss via absorption imaging of the atom 
cloud. The discrete number of trapped ions is clearly reflected in the dis- 
played histogram of atom numbers. 

number of ions in the trap. Up to five ions were trapped simultane- 
ously and detected with high fidelity. The atom loss rate increases 
nonlinearly with ion number mainly because the interionic repulsion 
prevents the ions from all occupying the trap center where the atomic 
density is maximal. While ion detection method 2 does not distinguish 
ionic masses, it has advantages in terms of experimental stability and 
does not require the trapping of ions amenable to laser cooling or other 
fluorescence based detection techniques. 



